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ABSTRACT 



The very young brown dwarf candidate Cha Ho- 8 was recently discovered to have a close (~1 AU) companion by means of radial 
velocity monitoring. We present additional radial velocity data obtained with UVES/VLT between 2007 and 2010, which significantly 
improve the orbit determination of the system. The combined data set spans ten years of radial velocity monitoring for Cha Ho- 8. A 
Kepler fit to the data yields an orbital period of 5.2 years, an eccentricity of e=0.59, and a radial velocity semi-amplitude of 2.4 km s"' . 
A companion mass M2 sin ; (which is a lower limit due to the unknown orbital inclination i) of 25±7 Mj^p and of 3 1±8 Mj^p is derived 
when using model-dependent mass estimates for the primary of 0.07 Mq and 0. lOM©, respectively. Assuming random orientation 
of orbits in space, we find a very high probability that the companion of Cha Ho- 8 is of substellar nature: With a greater than 50% 
, probability (/ > 60°), the companion mass is between 30 and 35 Mjup and the mass ratio M2/M1 smaller than 0.4; with a greater than 

. 87% probability (i > 30°) the companion mass is between 30 and 69 Mj^p and the mass ratio smaller than 0.7. The absence of any 

evidence of the companion in the cross-correlation function together with the size of the radial velocity amplitude also indicate a 
(— I ' mass ratio of at most 0.7, and likely smaller. Furthermore, the new data exclude the possibility that the companion has a mass in 

^ ^, the planetary regime (<13 Mjjip). We show that the companion contributes significantly to the total luminosity of the system: model- 

I . dependent estimates provide a minimum for the luminosity ratio L2/L1 of 0.2. Cha Ha 8 is the fourth known spectroscopic brown 

O ' dwarf or very low-mass stellar binary with determined orbital parameters, and the second known very young one. With an age of 

only about 3 Myr, it is of particular interest to very low-mass formation and evolution theories. In contrast to most other spectroscopic 
C/3 , binaries, it has a relatively long orbital period and it might be possible to determine the astrometric orbit of the primary and, thus, the 

■ orbital inclination. 
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^ 1. Introduction 2007); and 2MASS J03202839-0446358, an SBl at the M/L 

f~vj transition (Blake et al. 2008). Furthermore, for a few close vi- 

'-J Brown dwarf (BD) and very low-mass stellar (VLMS, M < ^^^j bd/VLM binaiies follow-up RV measurements permitted 

^ 0.1 Mo) bmanes ai-e key objects to understand formation additional consti-aints of their orbits, e.g., G7569B (Zapatero 

O and evolution in the low-mass regime, in particular when 0^^^^^ 2004; Simon, Bender & Prato 2006; Konopacky 

2 they are very young. Roughly a hundred BD/VLMS bina- 2OIO); and 2MASS J07464256+2000321 (Konopacky et 

7: ;ries have been detected m recent years, predommandy by al. 2010). When combined with angular distance measurements 

> , direct (adaptive optics or Hubble Space Telescope) imaging ^^j^p^g detections, spectroscopic binaries allow valuable dy- 

(e.g., |http://www. vlmbinaries.org). For several of these visual ^^^^^^^ ^^^^ determinations. Either the astrometiic orbit or 

BDA^LM binai-ies, it was possible to direcdy monitor the or- ^j^^ ^^j^p^g jigj^^ ^^^^^ ^^^^-^^^ inclination / of the or- 

^ bital motion and derive constraints on the orbital parameters and ^^^^^ pj^^^^ ^j^i^j^ j^i^^i^g parameter in the spectroscopic 

■ ■ mass (e.g., Bouy et al. 2004; Dupuy et al. 2009; Konopacky ^^^^ ^^^^^ • ^^^^ ^^^^ important input parame- 

et al. 2010; Stumpf et al. 2010). High-resolution spectroscopic evolutionary models, which at ages <10Myr and masses 

monitoring for radial velocity (RV) vanations led to the de- <o 3 rely on only three measurements, the two components 

tection of a few very close (<1 AU) specti'oscopic BDA'LM y^^^g gdipsing BD binary 2M0535-05 (Stassun et al. 

binaries. For four of them, RV orbital parameters were deter- 20O6; Mathieu et al. 2007) and the low-mass component of the 

mined: PPl 15, a double-hned spectroscopic binary (SB2) in the ain sequence star UZTauE (0.294±0.027Mo; Prato et al. 

Pleiades (Basri & Martin 1999); 2MASS J05352184-0546085 2002) 

(hereafter 2M0535-05), a very young eclipsing SB2 system in . 2M0^3S 05 ^haHryf^1 is the second known verv 

r>r;^„ (-Qtoco„„ \/rcth;»„ \/ci»„t; ru^u^R c ^"^"^ 2MU3^3-U3, t^naHQfjy is tne secona Known very 
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Orion (Stassun, Mathieu & Valenti 2006); Cha Ha 8, a very 
young spectroscopic binary in Chamaeleon I (Joergens & Miiller 



young (~3Myr) BD/VLM spectroscopic binary and, there- 
fore, of particular interest in the context of BD/VLM evo- 
lution and formation. Cha Ha 8 is a very low-mass member 



* Based on observations obtained at the Very Large Telescope of the 

European Southern Observatory at Paranal, Chile in program 279. C- 

5018(A), 380.C-06 19(A), 082.C-0023(A+B), 385.C-05 10(A). ' Simbad name: [NC98] Cha HA 8 
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Table 1. RV measurements of Cha Ho- 8. 



Date 


HJD 


RV 


ctrv 






[km s-'] 


[kms-'] 


2000 04 05 


2451639.61095 


14.591 


0.400 


2000 04 24 


2451658.72597 


15.177 


0.400 


2002 03 06 


2452339.68965 


17.499 " 


0.193 


2002 03 22 


2452355.65264 


17.355 " 


0.150 


2002 04 16 


2452380.61646 


17.651 " 


0.150 


2002 04 19 


2452383.57565 


17.451 " 


0.379 


2005 03 21 


2453450.62080 


14.791 " 


0.150 


2006 04 10 


2453835.65109 


16.035 


0.150 


2006 07 09 


2453926.50137 


16.403 


0.150 


2007 03 15 


2454174.66101 


17.358" 


0.402 


2007 03 22 


2454181.69756 


17.563 " 


0.333 


2007 06 08 


2454260.49431 


17.449 " 


0.182 


2007 07 21 


2454302.54665 


17.572" 


0.150 


2008 01 02 


2454467.82631 


17.573 " 


0.150 


2008 07 09 


2454657.48163 


18.190" 


0.150 


2009 01 02 


2454833.80683 


18.780" 


0.155 


2010 01 09 


2455205.74852 


13.900" 


0.150 


2010 03 01 


2455257.26535 


14.533 " 


0.203 



Notes. Listed are new RV data (bottom panel) and previous RV mea- 
surements (top panel), which have been re-processed for this work. HJD 
is given at the middle of the exposure; cr«i/ is the estimated error in the 
relative RVs. An additional error of about 400 m s"' has to be taken into 
account for the absolute RVs. 

RV value is the average of two (three in the case of 2010 03 01) 
single consecutive measurements. 

(M5.75-M6.5, Comeron et al. 2000; Luhman 2004, 2007) of 
the nearby (~160pc) Cha I star-forming region. An estimate 
of its mass based on evolutionary model tracks (BarafFe et al. 
1998) yields 0.07 Mq when using the values for effective tem- 
perature and bolometric luminosity by Comeron et al. (2000) 
and O.IOM0 when using those by Luhman (2007), respectively. 
Thus, Cha Ha 8 is either a BD or a VLMS. A more detailed de- 
scription of Cha Hot 8 and its properties can be found in Joergens 
& Muller (2007). 

RV data obtained until 2007 (Joergens & Muller 2007) 
showed that Cha Ha 8 has a companion in a few years orbit 
with a potentially very low mass close to or in the planetary 
mass regime. This made Cha Ha 8 a candidate for being the first 
BD/VLMS with an RV planet. We present here new RV mea- 
surements for Cha Ha 8 which significantly improve the RV or- 
bit for the system. 

2. Radial velocities and orbital solution 

Spectroscopic observations of Cha Ha 8 were carried out be- 
tween 2000 and 2010 with the Ultraviolet and Visual Echelle 
Spectrograph (UVES, Dekker et al. 2000) attached to the VLT 
8.2 m KUEYEN telescope at a spectral resolution AIM of 40 000 
in the red optical wavelength regime. RVs were measured from 
these spectra based on a cross-correlation technique employing 
telluric lines for the wavelength calibration. Details of the data 
analysis can be found in Joergens (2006; 2008). 

RV measurements for spectra taken between 2000 and 2002 
provided already evidence of an RV companion to Cha Ha 8 
(Joergens 2006). Follow-up RV monitoring between 2005 and 
March 2007 allowed the determination of a first spectroscopic 
orbital solution for Cha Ha 8 (Joergens & Muller 2007). Here, 
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2000 3000 4000 5000 

Time [JD-2450000] 

Fig, 1, RV measurements of Cha Ha 8 between 2000 and 2010 
based on UVES/VLT spectra. New RV data measured between 
Jun 2007 and Mar 2010 are marked with black-filled points, pre- 
viously obtained RV data with grey-filled points. Overplotted 
is the best-fit Keplerian orbit, which has a semi-amplitude of 
2.4km s ', a period of 5.2 years and an eccentricity of e=0.59. 



we present new RV data based on UVES spectra taken between 
June 2007 and March 2010. Table[T] lists the new RV measure- 
ments of Cha Ha 8 and previous RV data. The latter were re- 
processed for this work by using a slightly refined wavelength 
range for the cross-correlation to exclude a region that had been 
found to be contaminated by telluric lines and imperfectly cor- 
rected CCD cosmetic blemishes and that introduced additional 
scatter. The given errors ctrv are in most cases the sample stan- 
dard deviation of two consecutive individual measurements tak- 
ing into account an additional minimum error of 150 m s ^ This 
minimum error accounts for both RV noise due to activity and 
fluctuations in the standard deviation for a small number of mea- 
surements. It is derived from an estimate of the RV noise based 
on the night-to-night rms scatter in the data. 

The combined data set spans 9.9 years of RV monitoring for 
Cha Ha 8 between 2000 and 2010. We calculate an improved or- 
bital solution based on these data by fitting the six free parame- 
ters period, periastron time, eccentricity, longitude of periastron, 
RV semi-amplitude, and system velocity. The fitting procedure 
uses a global optimization method that is based on the genetic al- 
gorithm PIKAIA (Charbonneau 2002), as in the previous work. 
The reduced of the orbital fit is 1.45. 

The fitted Kepler orbit is that of a companion revolving 
Cha Ha 8 with a period of 1895d (5.2 yr) on an eccentric 
(e=0.59) orbit and causing an RV semi-amplitude of 2.4 km s"' . 
The derived lower limit to the semi-major axis of the primary 
a\ sin / is 0.34 AU. Figure[T]shows the RV measurements and the 
RV curve of the best-fit Kepler model. The complete list of de- 
termined orbital elements is given in Table|2l 

The minimum mass of the companion M2 sin / in the case of 
a single-lined RV orbit depends on the primary mass, which is 
not very precisely determined for Cha Ha 8, as is common in this 
mass and age regime. Using the two available estimates for the 
primary mass (0.07 and 0. IOM0, cf. Sect.lTJ, the mass M2 sin / of 
the companion is determined to be 25 and 31Mjup, respectively. 
The given errors in Mo sin / (cf. Table|2]l are based on the fit and 
do not take into account additional possible errors in the primary 
mass, such as those introduced by the evolutionary models used 
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to estimate the primary mass. The semi-major axis of the com- 
panion, which also depends on the primary mass estimate, is of 
the order of 1 AU. 

The new orbit solution for Cha Ha 8 is based on more epochs 
of RV measurements that cover a significantly larger portion of 
the orbital phase than the one derived by Joergens & Miiller 
(2007). As illustrated in Fig.[T] previous data covered only about 
half of the total orbit, while the new measurements monitor 
the system at orbital phases that were so far completely un- 
constrained. Thus, the constraints on the orbit are now much 
stronger than in 2007 (one should not be misled by the appar- 
ently superior goodness of the fit of the previous orbit, which 
is in fact a result of undersampling). Because of this much 
larger phase coverage, the refinement of the RV data, and the 
introduction of a minimum error, an orbital solution is derived 
that has a significantly longer orbital period, and a slightly 
larger RV semi-amplitude and eccentricity than previously an- 
ticipated. The newly derived companion mass is consequently 
higher but being within 1 cr still consistent with that found pre- 
viously (16+4 Mjup and 20+5 Mjup; note that the error ranges for 
M2 sin / are quoted incorrectly in Joergens & Miiller 2007). 

It has already been demonstrated in detail (Joergens & 
Miiller 2007) that the detected RV variations cannot be caused 
by activity. The main arguments are summarized in the follow- 
ing: (i) The timescale of the RV variations is much too long to 
be explained by any phenomena that causes rotational modula- 
tion (the rotation rate of Cha Ha 8 is of the order of a few days; 
Joergens & Guenther 2001; Joergens et al. 2003). (ii) Accretion 
can also be excluded as source of RV noise because neither mid- 
IR (Persi et al. 2000; Comerdn et al. 2000) nor L-band obser- 
vations (Jayawardhana et al. 2003; Luhman et al. 2008) indicate 
the presence of disk material surrounding Cha Ha 8. Ha mea- 
surements also provide no evidence of accretion (e.g. Mohanty 
et al. 2005). 



3. Discussion and conclusions 

Following up on the discovery of an RV companion orbit- 
ing the BD candidate Cha Ha 8 (Joergens & Miiller 2007), we 
have monitored the system between 2007 and 2010 with high- 
resolution spectroscopy with UVES at the VLT. The combined 
data set spans 9.9 years of RV monitoring of Cha Ha 8 and con- 
stitutes the longest RV monitoring program of a BD or VLMS. 
We have presented a significantly improved RV orbit solution 
of the system. The best-fit Kepler orbit of Cha Ha 8 has a 5.2 yr 
period, an RV semi-amplitude of 2.4kms \ and an eccentric- 
ity of 0.59. We determined the spectroscopic minimum mass 
Ml sin / for the companion to be 25+7 Mjup and 31+8Mjup for 
the two available model-dependent mass estimates for the pri- 
mary of 0.07 Mq and O.IOM0, respectively. 

Based on the new RV data and orbit solution, we have been 
able to exclude that the companion of Cha Ha 8 has a mass 
in the planetary mass regime (< 13Mjup). New RV measure- 
ments (black-filled points in Fig.[TJ significantly improve the 
phase coverage compared to previously available data (grey- 
filled points in Fig.[TJ and show that the RV amplitude is larger 
than previously anticipated and that the orbit is inconsistent with 
a planetary mass of the companion. 

The companion of Cha Ha 8 is with high probability of sub- 
stellar nature. This follows based on statistical considerations, 
as described in this paragraph, and it follows also independently 
based on the absence of companion lines, as described later 
While RV orbits permit only the determination of lower mass 



Table 2. Orbital and physical parameters derived for the best-fit 
Keplerian model of Cha Ha 8. 



Parameter 


Value 


/'(days) 


1895 ± 132 


T (HJD-2450000) 


3163 ± 214 


e 


0.59 ± 0.22 





106 + 26 


if (kms-') 


2.433 ± 0.44 


Vo(kms-') 


16.693 ± 0.75 


f{m) ilQ-' Me) . . 


1.488 


M2sin;(Mj„p) ... 


25 + 7°, 31 ± 8" 


ai sin i (AU) 


0.34 ± 0.09 


a2 (AU) 


1.02 ±0.06°, 1.17 ±0.07° 




18 


Span (days) 


3618 


a- (0-C)(m/s) .... 


228 




1.45 



Notes. The given parameters are: orbital period, periastron time, ec- 
centricity, longitude of periastron, RV semi-amplitude, system veloc- 
ity, mass function, lower limit of the companion mass, lower limit of 
the semi-major axis of the primary, semi-major axis of the companion, 
number of measurements, time span of the observations, residuals, re- 
duced X' ■ 

Derived parameter based on two available estimates for the primary 
mass of 0.07 Mq and 0. 10 Mq. No further errors of the primary mass, 
e.g. as introduced by the use of evolutionary models, have been taken 
into account here. 



limits M2 sin / due to the unknown inclination for orbits ori- 
ented randomly in space, there is a high probability that the true 
mass lies not too far from this limit. In the following, we calcu- 
late the absolute companion mass M2 for different inclinations. 
For this purpose, M2 is derived directly from the mass function 
rather than using the value of M2 sin / given in Table|2] which is 
for non-negligible companion masses valid only for large incli- 
nations. Furthermore, it is taken into account that the companion 
contributes increasingly to the total luminosity when considering 
increasingly smaller inclinations (see paragraph after next) and, 
therefore, that the primary mass becomes increasingly low com- 
pared to the value derived for the unresolved luminosity. Here, 
the luminosity determined by Luhman (2007) for the unresolved 
system is used, which corresponds to a mass of the unresolved 
source of 0. 10 Mq. We find that with a probability of greater than 
50% (inclination greater than 60°), the companion of Cha Ha 8 
has a mass between 30 and 35Mjup. With a probability higher 
than 87% (/ > 30°), it has a mass between 30 and 69 Mjup. 

The mass ratio q = M2/M1 of Cha Ha 8 seems to be on the 
small side compared to other known BD/VLM binaries: it is at 
minimum 0.3, with more than 50% probability smaller than 0.4, 
and with more than 87% probability smaller than 0.7. This is 
smaller than the mass ratio of the majority of BD/VLMS de- 
tected through direct imaging, where 68% have q>0.8 (value de- 
rived for binaries with a total mass <O.2M0 as listed in the VLM 
binary archive as of July 2009 and applying small updates; cf. 
also Burgasser et al. 2007). 

The companion of Cha Ha 8 contributes significantly to the 
total luminosity of the system. The luminosity ratio L2/L1 is es- 
timated for different inclinations by employing low-mass evo- 
lutionary model tracks (Baraffe et al. 1998) and assuming the 
coevality of the binary components. For the luminosity of the 
unresolved system and its effective temperature, the values de- 
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termined by Luhman (2007; Ltot=O.O37L0, reff=3024K) ai-e 
used, which correspond to a mass of the unresolved source 
of O.IOM0. At minimum, i.e. for /=90°, the luminosity ratio 
L2/L1 of the ChaHaS system is found to be 0.2. For larger 
inclinations, it becomes increasingly large, e.g., L2/Li=0.3 for 
;=60° and L2/L1-QJ for ;=30°, respectively. The limiting case 
is to consider a companion that is as luminous as the pri- 
mary (i.e., Li-L2-Liotl2). Comparison with evolutionary tracks 
shows that in this case the system would be slightly older 
(almost 5Myr) and both components would have a mass of 
Ml -M2 a;0.085 Mq. This constitutes the model-dependent mini- 
mum mass of the primary when using the parameters by Luhman 
(2007). Furthermore, it implies a minimum orbital inclination of 
24°. However, we show in the next paragraph that the compo- 
nents cannot have equal luminosities and masses but that L2lL\ 
and q are instead significantly smaller than 1 . 

There is no evidence of spectral lines of the companion in 
the spectra: neither exhibits the cross-correlation function a sec- 
ond peak, nor does its width or shape (bisector velocity span) 
correlate with RV. This implies that the companion is either 
too faint to be directly detected and/or that its signal is en- 
gulfed completely by the rotational broadening of the line profile 
(vsin/=16±3kms"'; Joergens & Guenther 2001). In the latter 
case, blending effects can cause the net RV signal to be weaker 
than the true RV amplitude. To determine the maximum mass 
and luminosity ratio that is consistent with the absence of sec- 
ondary lines and at the same time with the measured RV varia- 
tion, we simulate rotationally broadened Gaussian functions rep- 
resenting the two components of a hypothetical binary using the 
eccentricity and period as determined here, the primary mass of 
ChaHaS, allowing RV amplitudes of >2.4kms ', inclinations 
of 24-90°, and involving again evolutionary models to estimate 
the luminosity ratio. The case in which both components have 
equal masses, and, therefore, equal luminosities, can be excluded 
immediately because the net RV signal would then be zero. More 
generally, we find that the mass ratio cannot be too high because 
the simulated RV curves would be of too small net amplitude 
due to blending. We derive a conservative upper limit to botlQ 
the mass and luminosity ratio of about 0.7. 

Because of its very young age (~3 Myr) and the small num- 
ber of known BD/VLMS spectroscopic binaries with determined 
orbital parameters, the very low-mass binary system Cha Ho' 8 
is potentially of great interest for understanding formation and 
evolution close to and below the substellar borderline. Cha Ha 8 
has a relatively long orbit for a spectroscopic system and might 
therefore become one of the few spectroscopic binaries for 
which the relative astrometric orbit of the primary can be de- 
termined and for which absolute dynamical masses can be mea- 
sured. Therefore, an astrometric monitoring program of this ob- 
ject is ongoing. 
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^ Note that for young BDs, the mass-luminosity relation according 
to the applied evolutionary tracks is much less steep than for old field 
BDs, leading to constraints on the mass ratio that are of similar order as 
the luminosity ratio. 
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